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Nonextensive triplet in geological faults system 
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PACS 97 . 10 . Kc - Earthquakes 

PACS 97. lO.Yp - Star counts, distribution, and statistics 

PACS 05.90.+m- Other topics in statistical physics, thermodynamics, and nonUnear dynamical 

systems 

Abstract - The San Andreas fault (SAP) in the USA is one of the most investigated self- 
organizing systems in nature. In this paper, we studied some geophysical properties of the SAP 
system in order to analyze the behavior of earthquakes in the context of Tsallis's g^-Triplet. To 
that end, we considered 134,573 earthquake events in magnitude interval 2 < m < 8, taken 
from the Southern Earthquake Data Center (SCEDC, 1932 - 2012). The values obtained ("q- 
Triplet" ^{g'statj^senj^rei}) rcvcal that the ^stat-Gaussian behavior of the aforementioned data ex- 
hibit long-range temporal correlations. Moreover, qsen exhibits quasi-monofractal behavior with 
a Hurst exponent of 0.87. 
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Introduction. — Earthquakes are among the most 
complex spatiotemporal phenomena investigated in the 
context of self-organized criticality (SOC), introduced in 
Ref. [I]. In this regard, let us consider the so-called fault 
systems, a complex phenomenon related to the deforma- 
tion and sudden rupture of some parts of the Earth's crust 
driven by convective motion in the mantle. One of the first 
examples of self-organizing systems in nature 2 is the 
San Andreas Fault (SAF) in California. The SAP, one 
of the world's longest and most active geological faults, is 
~1200 Km long, ~15 Km deep, and about 20 million years 
old. It forms the boundary between the North American 
and Pacific plates and is classified as a right lateral strike- 
slip fault, although its movement also involves comparable 
amounts of reverse slip 3 . From the geophysical stand- 
point , a considerable number of investigations have been 
conducted in order to better understand the complexity of 
this system (see, e.g., |4 and references therein). In con- 
trast to the complexity of earthquakes, empirical laws are 
extremely simple, e.g. the Gutenberg-Richter law, which 
gives the number of earthquakes with a magnitude M > m 
[5] , and the Omori law for temporal distribution of after- 
shocks [gI. 

Several studies have demonstrated that seismicity ex- 



hibits an out-of-equilibrium behavior that is being inves- 
tigated by different authors, e.g. studies based on wavelet- 
based multifractal analysis 7 and nonextensive statistical 
mechanics [8 -10 , among others. In the present study, we 
consider nonextensive formalism, which is a generaliza- 
tion of Boltzmann-Gibbs statistical mechanics (B-G statis- 
tics) for out-of-thermal equilibrium systems and described 
by the entropic parameter q. The celebrated Boltzmann- 
Gibbs (B-G) statistics is recovered at g = 1 11 -13 . This 



parameter measures the degree of nonextensivity in the 
stochastic process. 

Tsallis statistics is based on the g'-exponential and q- 
logarithm, two central functions defined by 



exp r/) = [l + (l-g)/] 



1/1-q 
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which produces entropy Sq 13 , associated with q- 
statistics, 
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Fig. 1: Time series for magnitude of earthquakes along SAF. 
The peaks denote the maximum magnitudes. 



where the Boltzmann-Gibbs entropy, usual exponential 
and logarithm are recovered if g = 1. 

This theory has been successfully applied to many com- 
plex physical systems such as geological faults 10 and as- 



trophysical systems 14 16 . In 2004, Tsallis 17 proposed 
the existence of a three-parameter set {qstatAsenArei)i also 
known as (7- Triplet, characterized by metastable states in 
nonequilibrium, where qstat > 1, (Isen < 1 and Qrei > 1- 
When (fea^,gsen,^rez) = (1,1,1), the sct dcuotcs the B-G 
thermal equilibrium state. Burlaga and Vihas 18 used 
this triplet to describe the behavior of two sets of daily 
magnetic field strength performed by Voyager 1 in the so- 
lar wind in 1989 and 2002. In 2009, de Freitas and De 
Medeiros p6] presented a physical corroboration of the q- 
Triplet, based on analyses of the behavior of three sets 
of daily magnetic field strength observed by different so- 
lar indices. More recently, Ferri, Savio and Plastino 19 
showed a physical implication of this triplet for the ozone 
layer in Buenos Aires, Argentina. 

The main aim of this study is to analyze the behavior of 
physical parameters directly reflecting seismic activity in 
the context of Tsallis g'-Triplet's formalism, and to com- 
pare the properties of this g-Triplet with those expected 
for a metastable or quasi-stationary dynamical system de- 
scribed by nonextensive statistics. In this context, we fo- 
cus our attention on the magnitude-values for SAF data 
M{t) and their hourly variability dMT-{t). Following the 
ideas presented in Ref. |^, we focus our investigation on 
the "return" or fluctuation dM^(t) = M{t + r) - M(t), 
which denotes the differences between "avalanche" sizes 
obtained at time t -\- r and at time t. With respect to 
seismic activity, this analysis also checks the validity of 
g-Central Limit Theorem, the so-called g-CLT, recently 
conjectured by Umarov, Tsallis and Gell-Mann |21 . 

The remainder of this paper is organized as follows: in 
Section 2, we present our seismic sample; the main results 
and discussions are presented in section 3; and, finally, 
conclusions are put forth in the last section. 

The seismic data. — Fig. [T] shows the time se- 
ries for magnitude M of earthquakes along the SAF, in 
the interval 2 < M < 8, with 134,573 events. These 
were taken from the Southern California Earthquake Data 
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Fig. 2: Map of seismicity in SAF system, showing epicenters 
of earthquakes considered in this study (source SCEDC). 



Center (SCEDC) for 1932 to 2012. This range was cho- 
sen because for small magnitudes it has the limitation of 
seismic monitoring in the area, since many such events 
are unregistered. The lower panel of this figure shows 
non-overlapping magnitude fluctuations (return) in M for 
dMi {t) . Fig. [2] illustrates the distribution of events con- 
sidering the SAF map . 

Fig. [2] shows the data and the San Andreas fault sys- 
tem. This system is more than 800 miles long and extends 
to depths of at least 10 miles. The fault is a complex zone 
of crushed and broken rock ranging from a few hundred 
feet to a mile wide. Many smaller faults branch from and 
join the San Andreas fault zone. Almost any road cut in 
the zone shows a myriad of small fractures, fault gouge 
(pulverized rock), and a few solid pieces of rock 4 . The 
movement that occurs along the fault is a right-lateral 
strike-slip forming the tectonic boundary between the Pa- 
cific Plate and the North American Plate. 

Results and discussions. — In this sec- 
tion, we show results after the estimation of "g- 
Triplet"={gsta^, ^sen,<7rez} bascd on SAF data from 1932 
to 2012 (see Fig. fl]). These results are presented in three 
subsections, each associated to the properties of one of 
the g''s. 

On the behavior of the q-stationary parameter. For 
time series M{t)^ increment fluctuations due to its vari- 
ability over timescale r is given as dM^(t) = M{t + r) — 
M{t). The values of qstat are derived from probability 
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Fig. 3: Linear correlation between \nq[PDF] and [dMi(t)/cri]^, 
where Qstat = 1.364 ± 0.04, with R^ = 0.992 and x'^/dof = 
7.0236 X 10"^ 



distribution functions (PDFs). These PDFs are obtained 
from the variational problem using the continuous version 
for the nonextensive entropy given by eq. ([3| 



PDF = A,[l^{q-l)B,dMr{t)^Y 



(4) 



the entropic parameter q is related to the size of the tail in 
the distributions [15 
are given by 



and coefficients Aq and Bq for g > 1 
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for further details see Ref. i22i. 

Following the same procedure described by 19 , we var- 
ied the index q between 1.0 and 2.0, making a linear ad- 
justment in each computational iteration and evaluating 
the specific correlation coefficient R^ . The best linear 
fit is obtained for qstat = 1-364 ± 0.04 with R'^ = 0.992 
as shown in Fig. [3J It should be emphasized that this 
Qstat value is fully consistent with the bounds obtained 
from several independent studies involving the nonexten- 
sive Tsallis framework (see, e.g. [23^). The PDF for the 
return drM{t) on scale r = 1 is shown in Fig. [4j On 
this scale we can conduct a closer investigate of a pos- 
sible correlation between events . Our study used the 



Levenb erg-Mar quardt method 24 , 25 to compute PDFs 
with symmetric Tsallis distribution from Equation Q. In 



suits are consistent with the value expected for nonlinear 
systems, where the random variable is the sum of strongly 
correlated contributions 15,18 26. In this respect, we 



Fig. 4: In black circles is ploted distribution of increment for 
SAF data. In solid black line, the ^stat-Gaussian distribution 
based in eq. Q with Bq = 0.858 it 0.16. In dashed line, the 
best fit with a standard Gaussian. 

showed that PDFs for the return dMi{t) have fat tails 
with a g-Gaussian shape. 

On the behavior of the q-sensibility parameter. Values 
of the g^en-iiidex are directly related to system instability 
and entropy growth. These values can be obtained from 
multifractal (or singularity) spectrum /(a), where a is the 
singularity strength or Holder exponent. Spectrum f{a) 
is derived via a modified Legendre transform, through the 
application of the MF-DFA5 method [27]. This method 



consists of a multifractal characterization of a nonstation- 
ary time series, based on a generalization of detrended 
fluctuation analysis (DFA). MFDFA performs best when 
the signal is a noise-like time series. However, there is 
also difficulty in visualizing the difference between walk 

before 
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and noise-like time series. As suggested by 
application, it is necessary to run a DFA and verify if the 
value of the Hurst exponent is less than 1.2. For SAF 
data we obtain a Hurst exponent of 0.87, indicating that 
the MFDFA method can be employed directly without 
transformation of the time series. 

The (^sen-index denotes sensitivity at initial conditions. 
For present purposes, we used the expression defined by 
Lyra and Tsallis 29l for the relation between g^en and 
multifractality in dissipative systems, as follows: 



1 



1 -^s 



(7) 



where amin and amax denotes the roots of the best-fit. 

The multifractal characterization of these data is shown 
in Figure [s] These spectra /(a), calculated for SFA 
data, show a narrow Holder exponent interval with 



this adjustment, we found Bq = 0.858 ± 0.16. These re- a^ 



=0.924±0.04 and a^ 



=1.051±0.11. For multifrac- 



tal spectrum width, we obtained Aa=amax — Cimim re- 
sulting in a value of 0.127. Using Equation [7| we found 
that qsen = —6.647 ± 0.35. This negative value indicates 
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Fig. 5: The symbols are based on measurements of multifractal 
spectrum f{a) vs. a obtain from M{t). We obtain for SAF 
amin = 0.924±0.04 and amax = 1.051±0.11 a qsen = -6.747± 
0.35. The curve represent the best ajustment using a cubic fit 
to the data. 



that its distribution exhibits weak chaos [17 in the full 



dynamical space of the system 17,18 . Furthermore, this 



figure reveals that the behavior of our sample is similar to 
that of a monofractal-like time series. 

On the behavior of the q-relaxation parameter. The 
value of Qreh which describes a relaxation process, can be 
computed from an autocorrelation coefficient as a function 
of scale r defined by 



C{t) 



{[Siu + T)-{SiU))][SiU)-{SiU))]) 

{[S{u) - {S{uw) 



(8) 



In agreement with Tsallis statistics, we can estimate the 
value of Qrei by best fit on In^ C{r) vs. scale r, as shown in 
Fig. [g] (upper panel), where C{r) is given by Equation ([8|. 
In the nonextensive theory, this coefficient should decay 
following a power law, with increasing r, where slope s is 
given by s = l/{l — qrei)- From this adjustment, we obtain 
Qrei = 2.69 ±0.13 for SAF data. Moyano 30 suggests that 
the above procedure for calculating qrei only be used to 
describe stochastic processes with linear correlations. In 
other words, autocorrelation coefficient C{r) is not a good 
alternative to conveniently describe the non-linearity of a 
sample 16 . 

On the other hand, in B-G statistics, in contrast to 
the nonextensive theory, coefficient C{r) should decrease 
exponentially with an increasing r, following a C(r) = 
Aiexp(— r/ti) + ^2 exp(— r/t2) relation, with ti and ^2 
corresponding to the correlation or relaxation times. The 
fit shown in Fig. [6] (lower panel) reveals that ^2 ^ ti. As 
mentioned by 22 , this behavior is related to local equi- 
librium, and then a much slower decay for larger r. In 
agreement with these authors, this constitutes a neces- 
sary condition for the application of the superstatistical 
model, as described in Ref. [31j . 



Fig. 6: Upper panel: lug of autocorrelation coefficient C{r) vs. 
time delay r for SAF data. Lower panel: symbols represent 
the autocorrelation function for our sample and the gray line 
represents a double exponential fit with characteristic times of 
ti = 8.42 and t2 = 313.74 yielding a ratio about 37 between 
these two time scales (R^ = 0.964, x^/dof = 1.4 x 10~^ and 
time is expressed in order of hours). 



See 32 for further details and an extensive discussion 



about the estimation of Tsallis g-triplet. 

Conclusions. — We used a new approach to nonex- 
tensive formalism for hourly measurements of earthquakes 
along the SAF from 1932 to 2012. From these data 
we were able to estimate the values of the nonexten- 
sive three-index. We found that qstat = 1.364 ± 0.04, 
Qsen — —6.647 ± 0.35 and qrei = 2.69 ± 0.13. It is impor- 
tant to underscore that the result of the Qstat is consistent 
with the upper limit q < 2 obtained from several inde- 
pendent investigations |23J. In addition, the values of this 
triplet confirm the general scheme qsen < 1 < Qstat < Qrei , 
according to the nonextensive scenario proposed by Tsal- 
lis 
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These results reveal that this system is consis- 
tent with a nonequilibrium state, strongly suggesting that 
long-range correlations exist among the random variables 
involved in the physical process that controls seismic ac- 
tivity. 

Finally, it is worth mentioning that the nonextensive 
three-index can be recalculated by considering a spa- 
tiotemporal analysis for earthquakes along the SAF. This 
issue will be addressed in a forthcoming communication. 



Research activity at the Stellar Board of the Federal 
University of Rio Grande do Norte (UFRN) and Federal 
Institute of Rio Grande do Norte (IFRN) are supported 
by continuous grants from CNPq and FAFERN Brazilian 
agency. The authors would like to thank Cesar Garcia 
Pavao for their many helpful with maps this work. 



p-4 



Nonextensive triplet in geological faults system 



REFERENCES 

[I] Bak p., Tang C. and Wiesenfeld K., Phys. Rev. Lett., 
59 (1987) 381. 

[2] RuNDLE J. B., TiAMPO K. P., Klein W. and Martins 
J. S. S., Proc. Natl. Acad. Sci. USA, 99 (2002) 2514. 

[3] Wallace R. E. (Editor), The San Andreas Fault Sys- 
tem, California (United States Government Printing Of- 
fice, Washington) 1990. 

[4] SCHULZ S. S. and Wallace R. E., The San Andreas 
Faults, Jittp://puhs. usgs.gov/gip/earthq3/endnotes. html\ 
(1997) . 

[5] Gutenberg B. and Richter C. P., Bull. Seismol. Soc. 
Am., 34 (1944) 185. 

[6] Omori p., J. Coll. Sci. Imp. Univ. Tokyo, 7 (1894) 111. 

[7] Enescu B., Ito K. and Struzik Z. R., Annals of Disas. 
Prev. Res. Inst, 47 (2004) 16. 

[8] SoTOLONGO-CosTA O. and Posadas A., Phys. Rev. Lett, 
92 (2004) 048501. 

[9] SiLVA R., Pranca G. S., Vilar C. S. and Alcaniz J. S., 
Phys. Rev. E, 73 (2006) 026102. 

[10] Vilar C. S., Pranca G. S., Silva R. and Alcaniz J. 
S., Physica A, 377 (2007) 285 

[II] TsALLis C., J. Stat Phys., 52 (1998) 479. 

[12] Abe S. and Okamoto A. (Editors), Nonextensive Statis- 
tical Mechanics and Its Applications (Springer-Verlag, Hei- 
delberg) 2001. 

[13] Gell-Mann M. and Tsallis C. (Editors), Nonexten- 
sive Entropy- Interdisciplinary Applications (Oxford Univ. 
Press, New York) 2004. 

[14] BuRLAGA L. P. and ViNAS A. P. J., Ceophys. Res., 109 
(2004) 12107. 

[15] BuRLAGA L. P., Ness N. P. and Acuna M. H., ApJ, 
691 (2009) 82. 

[16] DE Preitas D. B. and De Medeiros J. R., Europhys. 
Lett, 88 (2009) 19001. 

[17] Tsallis C., J. Stat Phys., 52 (2004) 479 

[18] BuRLAGA L. P. and ViNAS A. P., Physica A, 356 (2005) 
375. 

[19] Perri G. L., Reynoso Savio M. P. and Plastino A., 
Physica A, 389 (2010) 1829. 

[20] Caruso P., Pluchino A., Latora V., Vinciguerra S. 
and Rapisarda A., Phys. Rev. E, 75 (2007) 5101. 

[21] Umarov S., Tsallis C. and Steinberg S., Milan J. 
Math., 76 (2008) 307. 

[22] QuEiROS S. M. D., Moyano L. G., de Souza J. and 
Tsallis C., Eur. Phys. J. B, 55 (2007) 161. 

[23] BOGHOSIAN B. M., Braz. Journ. Phys., 29 (1999) 91; 
Hansen S. H., Egli D., Hollenstein L. and Salzmann 
C., New Astronomy, 10 (2005) 379; Silva R., Alcaniz J. 
S. and Lima J. A. S., Physica A, 356 (2005) 509; Liu B. 
and GOREE J., Phys. Rev. Lett, 100 (2008) 055003; Car- 
VALHO J. C., Jr. do Nascimento J. D., Silva R., De 
Medeiros J. R., Astrophys. J. Lett, 696 (2009) 48. 

[24] Levenberg K., Q. Appl. Math., 2 (1944) 164. 

[25] Marquardt D., SI AM J. Appl. Math., 11 (1963) 431. 

[26] TiRNAKLi U., Beck C. and Tsallis C., Phys. Rev. E., 
75 (2007) 106. 

[27] Kantelhardt J. W., Zschiegner S. A., Koscielny- 
BuNDE E., Havlin S., Bunde a. and Stanley H. E., 
Physica A, 316 (2002) 87. 

[28] Eke A., Hermann P., Kocsis, L., Kozak L. R., Phys- 



iol. Meas., 23 (2002) 1. 
[29] Lyra M. L. and Tsallis C., PRL, 80 (1998) 53. 
[30] Moyano L., Thesis, CBPF, (1944) . 
[31] Beck, C. and Cohen, E. G. D., Physica A, 322 (2003) 

267. 
[32] Pavlos G. P., Karakatsanis L. P., Xenakis M. N., 

Sarafopoulos D. and Pavlos E. G., Physica A, 391 

(2012) 3069. 



p-5 



